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Abstract We collect the 2LAC and MOJAVE quasi-simultaneous data to investigate the 
radio-7 connection of blazars. The cross sample contains 166 sources. The statistic anal- 
ysis based on this sample confirms positive correlations between these two bands, but the 
correlations become weaker as the 7-ray energy increases. The statistic results between 
various parameters show negative correlations of 7-ray photon spectral index with 7-ray 
loudness for both FSRQs and BL Lacertae objects, positive correlations of 7-ray variabil- 
ity index with the 7-ray loudness for FSRQs, a negative correlation of the 7-ray variability 
index with the 7-ray photon spectral index for FSRQs, and negative correlations of 7-ray 
photon spectral index with 7-ray luminosity for FSRQs. These results suggest that the 
7-ray variability may be due to changes inside the 7-ray emission region like the injected 
power, rather than changes in the photon density of the external radiation fields, and the 
variability amplitude tends to be larger as the 7-rays are closer to the high energy peak 
of spectral energy distribution. No correlation of variability index found for BL Lacertae 
objects implies that variability behavior may differ below and above the peak energy. 

Key words: BL Lacertae objects: general - galaxies: jets - quasars: general - radio con- 
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1 INTRODUCTION 

Blazars, the most extreme class of Active Galactic Nuclei (AGNs), are observed at almost the full elec- 
tromagnetic spectrum from radio to 7-ray band. It contains two subclasses called BL Lacertae objects 
(BL Lacs) and Flat Spectrum Radio Quasars (FSRQs). The radiation from blazars is believed to be non- 
thermal emission of relativistic jets with small view angles to the line of sight (Urry & Padovani |1995t . 
The broadband spectral energy distribution (SED) shows two components in the v — vF v diagram. The 
lower energy component that peaks in infrared to X-ray band is believed to be produced by synchrotron 
process of relativistic electrons. The radiation mechanism of higher energy component is more debat- 
able until now. It is generally believed that the 7-ray emission is produced by the inverse Compton (IC) 
process of the same electron population that emits the synchrotron emission (e.g., Maraschi, Ghisellini 
& Cellotti 1992] Sikora, Begelman & Rees 1994). This indicates a connection between radio and 7-ray 
emission, but the connection could be indirect because the emission regions of radio and 7-ray may 
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be different due to synchrotron self-absorption. Moreover, some observations found the connection be- 
tween these two bands, at least during the flare period, e.g., the 7-ray flares occur after the ejections of 
new superluminal components (Jorstad et al. 20011. However, the location of 7-ray emission region and 
the source of seed photons are still uncertain (Marscher et al. 120 101 Leon-Tavares et al. 1201 11 Tavecchio 
et al. [20101 , and why only a small fraction of blazars are 7-ray loud is still an open question (e.g., 
Kovalev et ai. |2009l Linford et al |20Tfl Sikora et al. 120021 . 

The connection between radio and 7-ray emission can help us constrain the radiation process and 
the emission region in the jets, especially for 7-ray band. This has been considered since the EGRET 
era, but no confirmed conclusion had been made due to the limits of the sensitivity of the telescope 
and the erratic sample (Kovalev et al. |2009i >. Thanks to the high sensitivity, broad energy range, and 
large view field of Large Area Telescope (LAT, Atwood et al. 2009) on board the Fermi Gamma Ray 
Space Telescope (Fermi) launched successfully in 2008, much larger sample and more accurate photon 
flux and spectrum can be achieved. Many papers investigated the connection of these two bands, and 
almost all of them confirmed a positive correlation in Fermi blazars (Kovalev et al. 120091 Ghirlanda et 
al. 120101 1201 11 Linford et al 1201 II 120121 Ackermann et al. 1201 1 at . However, these previous results 
showed many different behaviors in different subclasses or energy bands. Moreover, these studies only 
used non-quasi-simultaneous data, or only studied the correlation of apparent fluxes which depend on 
distance. 

After two years survey of LAT, the second LAT AGN catalog (2LAC) gives a more accurate clas- 
sification and association of 1121 AGNs in 7-rays, and 886 sources in its clean sample (Ackermann 
et al. 1201 lb) . Some radio monitor programs are processing in Fermi era at different frequency, e.g., 
MOJAVE (Monitor of Jets in AGN with VLBA Equipment, Lister et al. |2QQ91 > program at 15GHz, 
the OVRO (Owens Valley Radio Observation, Richards et al. 120111 1 program at 15GHz, UMRAO (the 
University of Michigan Radio Variability Program) at 4.8, 8.0, and 14.5GHz, and so on. These pro- 
grams give us a good chance for investigating the connection between radio and 7-ray emission with 
quasi-simultaneous data. 

Abdo et al. ( 2012 ) obtained the 7-ray photon spectral index fitted with a single power-law and the 7- 
ray variability index, which shows the level of 7-ray variability. The spectral index is often discussed in 
many researches, as it reflects the peak of spectral energy distribution (SED, e.g., Ghisellini, Maraschi 
& Tavecchio 120091 Lister et al. 1201 U . It is possible that the variability index is also important for 
researching the 7-ray emission region and the origin of 7-ray variability. Ackermann et al. (1201 la) has 
discussed it in 7-ray range for 2LAC clean sample. 

In this paper, we statistically analyze both flux and luminosity correlations with the quasi- 
simultaneous data of 2LAC and MOJAVE for the first time. The data of seven different 7-ray energy 
bands are used to study the dependence of the radio-7 connection on 7-ray energy bands. The 7-radio 
luminosity ratio is defined as a parameter of 7-ray loudness to discuss the 7-ray spectrum and variability 
in detail (Arshakian et al. 120121 Lister et al. 1201 1) . Then we examine the correlations of 7-ray photon 
spectral index and 7-ray variability index with 7-radio luminosity ratio and other parameters in our 
sample. 

In section 2, we describe our sample and data, section 3 is the methods we use and the results, 
section 4 is discussion and summary. In this paper, we use a ACDM cosmology model with h=0.71, 
£l m =0.21, S1a=0.73 (Komatsu et al. 120091 ). We define the radio spectral index as S(v) oc v~ a , and the 
7-ray photon spectral index as dN/dE oc E~ r . 

2 SAMPLE AND DATA 

We choose the sources of 2LAC clean sample which have been monitored by the MOJAVE program. 
For the quasi-simultaneous data, we require a time interval matching the Fermi data (from 2008-08-04 
to 2010-08-01, Abdo et al. 120 121 ). Considering the gap of the radio monitor, and the probable variability 
time lag in these two bands, we use a broad time interval for radio data during 2008-01 to 201 1-04. In 
order to perform the K-correction and luminosity calculations, we reject the sources without redshift 
data. Finally, our sample contains 166 sources, with 127 FSRQs and 39 BL Lacs (hereafter cross sam- 
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pie). For FSRQs, there are 1 19 LSPs, and 8 sources without SED classification. For BL Lacs, there are 7 
HSPs, 13 ISPs, 18 LSPs, and 1 source without SED classification. The classification and redshift data of 
our cross sample are from 2LAC (Ackermann et al. 1201 lb) . The K-corrections for radio flux and 7-ray 
photon flux are calculated as 

s U0 =s' V0 {i + zY +1 (l) 

and 

N(E U E 2 ) = N'{E 1 ,E 2 )(l + z) r -\ (2) 
respectively. The radio and 7-ray luminosity are calculated by the observed flux as 

L Vo = 4 x 10- 26 ird 2 L Sl (1 + z)"- 1 erg s" 1 Hz- 1 (3) 
and (Ghisellini, Maraschi & Tavecchio |20091 > 

L(y u v 2 ) = Aird 2 L S\ VuV2 ){l + z) T - 2 erg s~\ (4) 

respectively. In equations (3) and (4), is the luminosity distance, and S v 2 ) is the 7-ray energy 
flux observed between frequencies v\ and v 2 which correspond to the energies E\ and E 2 . S (v\ ,is 2 ) is 
derived as 

S(y 1 ,v 2 ) = (2_r)[(^)i-r-i] erg s cm ^ ( ) 

Due to the absence of quasi-simultaneous data of radio spectral index, and generally the flat radio spec- 
trum for blazars, we assume that the radio spectral index is equal to for the sources in our sample. 

For Fermi data, we use the lGeV-lOOGeV photon flux in 2LAC (Ackermann et al. 1201 lbl l. and the 
100MeV-300MeV, 300MeV-lGeV, lGeV-3GeV, 3GeV-10GeV data i n 2FGL (Abdo et al. l20T2b . For 
upper limits, we use them for our analysis directly (Ackermann et al. 12011 al l, but we remove lOGeV- 
1 OOGe V data because too many sources only get upper limits in this band. We obtain the 1 OOMe V- 1 Ge V 
photon flux combined with the 100MeV-300MeV and 300MeV-lGeV data, the lOOMeV-lOOGeV pho- 
ton flux combined with the 100MeV-300MeV, 300MeV-lGeV and lGeV-lOOGeV data. For MOJAVE 
data, we use the mean value of all the VLBA observations for each source. 



3 RESULT 

We use the nonparametric Spearman test for our correlation analysis except that the luminosity corre- 
lation is tested by partial correlation analysis to remove the dependence on redshift. In the paper, we 
define a significant correlation when the chance probability is less than 0.001. 

The data of 15GHz radio flux and the 7-ray photon flux (lOOMeV-lOOGeV) are plotted in figure 1. 
The Spearman test confirms the positive correlation of radio and 7-ray flux, though the correlation has 
certain dispersion. This means that the sources with higher radio flux always have higher 7-ray flux. The 
correlation for BL Lacs seems to be better than that for FSRQs. The correlation coefficient and chance 
probability are 0.38 and 1.1 x 10~ 5 for FSRQs, and 0.44 and 0.006 for BL Lacs, respectively. 

Then we examine the correlation between the photon flux of individual 7-ray energy band and the 
radio flux. The correlations become worse when the energy band increases, especially for BL Lacs. This 
indicates that the connection between radio and 7-ray emission weakens as the 7-ray energy increases. 
There are no correlations of the radio flux with those of lGeV-3GeV, 3GeV-10GeV, and lGeV-lOOGeV 
for BL Lacs. Otherwise, for other bands including lOOMeV-lOOGeV, the correlation coefficients of BL 
Lacs are larger than those of FSRQs. For the radio and gamma-ray luminosity, we find positive correla- 
tions in all energy bands with the correlation coefficients of BL Lacs larger than those of FSRQs (Figure 
2). Correlation coefficients decrease when 7-ray energy gets higher. Figure 3 shows the changes of cor- 
relations with four energy bands, 100MeV-300MeV, 300MeV-lGeV, lGeV-3GeV and 3GeV-10GeV, 
respectively. 
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Fig. 1 15GHz radio flux versus 7-ray pho- 
ton flux (lOOMeV-lOOGeV). The filled circles 
represent FSRQs and the open ones represent 
BL Lacs. 



Fig. 2 15GHz radio luminosity versus 7-ray 
integral luminosity (lOOMeV-lOOGeV). The 
symbols are same as figure 1 . 
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Fig. 3 Changes of correlations with energy bands. The black filled squares represent cor- 
relation coefficients for FSRQs, the blue filled circles represent correlation coefficients for 
BL Lacs, the red filled triangles represent correlation coefficients for all sources in our cross 
sample, while the open ones represent the chance probabilities, respectively. Left shows the 
change of flux correlations, and the dot-dashed line represents 0.001 of the chance probabil- 
ity. Right shows the change of luminosity partial correlation coefficients. The corresponding 
chance probabilities are almost all less than 0.001, only one is 0.001. 



Figure 4 shows the plot of the 7-ray photon spectral index versus the 7-radio luminosity ratio 
(lOOMeV-lOOGeV). BL Lacs and FSRQs have a similar range in 7-ray loudness, and both of them 
have negative correlations between the 7-ray photon spectral index and the 7-ray loudness in most of 
7-ray energy bands, with the correlations for BL Lacs better than those for FSRQs. For FSRQs, we find 
positive correlations between the 7-radio luminosity ratio and 7-ray variability index (Figure 5). The 
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7-ray variability index is negatively correlated with the 7-ray photon spectral index for FSRQs (Figure 
6). 
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Fig. 4 The 7-radio luminosity ratio (7-ray lu- 
minosity is from lOOMeV to lOOGeV) versus 
7-ray photon spectral index. The symbols are 
same as figure 1 . 



Fig. 5 The 7-radio luminosity ratio (7-ray lu- 
minosity is from lOOMeV to lOOGeV) versus 
7-ray variability index. The symbols are same 
as figure 1. 
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Fig. 6 The 7-ray photon spectral index versus 
7-ray variability index. The symbols are same 
as figure 1. 



Fig. 7 The 7-ray integral luminosity 
(lOOMeV-lOOGeV) versus 7-ray photon 
spectral index. The symbols are same as 
figure 1. 



We also examine the correlations of radio flux and luminosity versus 7-ray photon spectral index 
and 7-ray variability index, respectively. We find the same trends for flux and luminosity. They have 
good correlations between radio flux and luminosity versus 7-ray photon spectral index for BL Lacs, 
respectively. There are no correlations for FSRQs. There are no correlations between the radio flux and 
luminosity and the 7-ray variability index (see Table 1). For the 7-ray luminosity of individual energy 
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band versus the 7-ray photon spectral index and the 7-ray variability index, we find positive correlations 
between the 7-ray luminosity and the variability index for FSRQs (see Table 1), faint positive corre- 
lations between the 7-ray luminosity and the photon spectral index for BL Lacs, and weak negative 
correlations between the 7-ray luminosity and the photon spectral index for FSRQs, especially in GeV 
band (see Figure 7). 

All the correlation coefficients and chance probabilities are listed in table 1 except for the case of 
the radio flux versus the 7-ray photon spectral index and the 7-ray variability index. 

4 DISCUSSION 

For our cross sample, we make the Kolmogorov-Smirnov (K-S) tests of the redshift distribution com- 
pared with the 2LAC clean sample and a full blazar sample, i.e., the bzcat blazar sample (Massaro et 
al. |20091 l. There is an intrinsic difference of the redshift distribution between the cross sample and the 
bzcat sample (D = 0.184, probability = 0.0005 for FSRQs, and D = 0.291, probability = 0.003 for BL 
Lacs). There is no obvious difference between 2LAC clean sample and our cross sample (D = 0.045, 
probability = 0.991 for FSRQs, and D = 0.159, probability = 0.366 for BL Lacs, see Figure 8). This 
means that there is an intrinsic difference in redshift distribution between 2LAC clean sample and bzcat 
sample, which may be caused by the limited number of 7-ray detected sources, or the special physics 
mechanisms for 7-ray sources. 




Fig. 8 The reshift distributions of three samples. Left is a full blazar sample as bzcat. Center 
is 2LAC clean sample. Right is our cross sample. 

There is no difference for FSRQs in 7-ray photon spectral index distributions of our cross sample 
and 2LAC clean sample (D = 0.068, probability = 0.78). For BL Lacs, the sources in our sample tend 
to have soft 7-ray spectrum than those in 2LAC clean sample (see Figure 9). This may be caused by a 
small fraction of HSP BL Lacs in our sample, because HSPs are generally weaker than ISPs and LSPs 
at the radio band (Ackermann et al. 1201 lal . For all the sources of our sample, they tend to have higher 
variability index than the 2LAC sources (Figure 10). This means that our sample selected the more 
variable sources and the more bright sources in 7-ray band because the variable sources determined by 
variability index must be both intrinsically variable and sufficiently bright (Ackermann et al. 1201 lbb . 

Ackermann et al. (1201 lab found the different behaviors of different energy bands for flux correla- 
tion. They found that BL Lacs exhibited larger p values than FSRQs in the full sample, but a opposite 
trend in OVRO sample. In OVRO sample, no correlation were found for BL Lacs above lGeV, which 
did not appear in the full sample. The different results of these three samples may be caused by the se- 
lection effects such as different redshift distributions or different proportions of HSPs, ISPs, and LSPs in 
different samples. So we also remind that the radio-7 connection, especially the flux correlation, should 
be used very carefully. 

The luminosity correlations between radio and individual 7-ray energy band are confirmed by our 
sample. However, both flux and luminosity correlation coefficients decrease as the 7-ray energy in- 
creases. For HSPs, the 15GHz and the GeV emission would be produced by low energy electrons, and 
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Fig. 9 The 7-ray photon spectral index distributions of two samples. Left is 2LAC clean 
sample. Right is our cross sample. 




Fig. 10 The 7-ray variability index distributions of two samples. Left is 2LAC clean sample. 
Right is our cross sample. 



there would be good correlations in these bands. But for ISPs and LSPs, the emission of low energy elec- 
trons could only extend to MeV band, and the correlations of radio and GeV bands would be weaker 
than those of radio and MeV bands. This would be the main reason why the correlation coefficients 
decrease as the 7-ray energy increases for both BL Lacs and FSRQs. Moreover, other reasons may also 
generate these trends, such as, other emission contents or emission regions in the GeV band. 

Lister et al. ( 1201 11 1 found a correlation between 7-ray loudness and 7-ray photon spectral index 
for BL Lacs, but no correlation for FSRQs. In our results, except for the energy below lGeV, both BL 
Lacs and FSRQs have negative correlations between the 7-ray loudness and the 7-ray photon spectral 
index. According to the average SEDs of blazar sequence (Fossati et al. 1998]), as the 7-ray spectrum 
becomes harder, i.e., the 7-ray photon spectral index decreases, the 7-ray peak of SED shifts to the 
higher energy. At the same time, in the energy band of Fermi/LAT, the ratio of the 7-ray luminosity to the 
radio luminosity would increase significantly. Then it is expected that the 7-ray loudness is negatively 
correlated with the 7-ray photon spectral index. This expectation is consistent with our results. This trend 
is more remarkable among different subclasses for BL Lacs. For FSRQs, it becomes more remarkable 
as the energy is larger than lGeV (see Table 1). This is due to the effective cooling above the SED peak. 

The correlations between the 7-ray variability and the 7-radio luminosity ratio suggest that the 
variability, at least the main variability may originate inside the 7-ray emission region like the changes 
of the injected power, because the changes in photon density of external radiation field such as variability 
of accretion disk could not affect the 7-ray loudness. This is opposite to the view of Paggi et al. ( 1201 \\ . 
These different trends of BL Lacs and FSRQs imply that variability behavior may distinguish below and 
above the peak energy (Ulrich, Maraschi & Urry 11997] Ackermann et al. 1201 lbt . As the 7-ray photon 
spectral index decreases, the high energy component shifts to the higher energy, and the 7-rays in Fermi 
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energy range becomes closer to the high energy peak of SED for FSRQs. So the negative correlation of 
the 7-ray photon spectral index with the 7-ray variability index for FSRQs indicates that the variability 
amplitude tends to be larger as the 7-rays are closer to the high energy peak of SED. 

The correlations of radio and 7-ray luminosity with 7-ray photon spectral index for BL Lacs reflect 
the blazar sequence among different subclass of BL Lacs (Fossati et al.[T998 ). For FSRQs, the negative 
correlation of 7-ray luminosity with 7-ray photon spectral index weas also found in Ackermann et al. 
( 1201 lbl) . where the correlation was very weak. In our results, the correlations become significant when 
the energy is larger than lGeV. These correlations would emerge if the more powerful sources have 
harder electron spectrums. However, the fast falling above the high energy peak of SED could make the 
integral luminosity of 7-ray smaller for the sources with the steeper 7-ray spectrum, even that the peak 
luminosity or the total luminosity is higher. This is consistent with the trend that the higher energy bands 
have better correlations (see Table 1). Our results, which seem opposite to the expectation of the blazar 
sequence, are caused by this reason with high possibility. So as mentioned in Chen & Bai ( 1201 II) . the 
7-ray luminosity in Fermi range can not be used simply to investigate the blazar sequence instead of the 
the peak luminosity or the total luminosity, at least for FSRQs. The same trend regarded as harder when 
brighter is always observed in single sources, e.g., for 3C454.3 (Abdo et al. 1201 II ). 

In summary, the radio and 7-ray emission are in good connection in Fermi/LAT blazars, but these 
correlations become worse as the 7-ray energy increases. Moreover, the flux correlations would be 
affected strongly by the selection effects. The 7-ray variability index is correlated with other parameters 
for FSRQs. These correlations suggest that the 7-ray variability may be due to changes inside the 7- 
ray emission region like the injected power, rather than changes in the photon density of the external 
radiation field, and the variability amplitude tends to be larger as the 7-rays are closer to the high 
energy peak of SED. The different variability behaviors below and above the peak energy may cause 
the different trends of BL Lacs and FSRQs. The negative correlations of the 7-ray luminosity with 7- 
ray photon spectral index suggest that the 7-ray luminosity in Fermi range can not be used simply to 
investigate the blazar sequence instead of the the peak luminosity or the total luminosity, at least for 
FSRQs. However, our results are also limited by the source number of our sample, especially for BL 
Lacs. These results will be tested with larger samples. 
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Table 1 The Results of Correlation Test, p is the correlation coefficient of Spearman test, r is o 
the correlation coefficient of partial correlation analysis, and P is the chance probability. 
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